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Abstract 
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> 

The possible cosmological effects of primordial fluctuation corrections 
. to the evolution equation of matter obtained from the Wheeler-De Witt 

' equation are explored. In particular, both the metric and a scalar matter 

field are expanded around their homogeneous values and the corrections 
induced on the scalar field fluctuation spectrum are perturbatively esti- 
mated. Finally, results of a preliminary numerical simulation to investi- 
O ■ gate the effects on large-scale structure formation are presented. 

cr 

I ■ The matter-gravity system may be studied quantum-mechanically through 

the canonical quantization of gravity within the superspace approach Q Q. 
One then obtains the Wheeler-De Witt (WD) equation and a matter-gravity 
wave function which, in analogy with the Born-Oppenheimer (BO) studies of 
■ molecules ||, may be factorized into two parts: one involving only gravitational 

degrees of freedom and the other involving both the gravitational and the matter 
degrees of freedom. Correspondingly, the WD equation in which time is absent 
is split into two pieces: one describing gravitation in an effective potential given 
by the mean energy-momentum tensor of matter and the other describing the 
matter whose evolution is parametrized by time which is derived from the semi- 
classical approximation to the gravitational wave function M |j| . The above is 
contingent on the hypothesis that the Planck mass is much larger than the mass 
of any matter field or any inverse length scale used to describe matter. 

The above approach has been examined in detail within the context of a 
minisuperspace model with matter j(| with the view of understanding under 
what conditions the BO approximation is valid and quantum cosmology leads 
to the usual physics (Schrodinger equation for matter and Einstein classical 
equations for gravity on scales larger than the Planckian) . It was found that such 
was the case in an inflationary scenario S and after ten or more Planck times the 
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usual physics ensues. Such a result was obtained by estimating perturbatively 
the fluctuation corrections due to quantum gravitational effects both on the 
equation of motion for the purely gravitational part and that for matter, since 
the actual solution of the coupled nonlinear equations is extremely complicated. 

The fluctuation corrections considered were a consequence of matter-gravity 
forming a closed system: fluctuations of gravitational origin generate corre- 
sponding fluctuations in the evolution of matter and it is such fluctuations that 
are neglected in the BO approximation. Henceforth we shall refer to them as 
fluctuation corrections. 

It has been previously suggested that inflation could in principle provide a 
causal mechanism for the origin of density perturbations that later grow to form 
large-scale structures 0, thus we feel it is of interest to explore the influence 
of the above fluctuation corrections in such a context. In order to do so we 
generalize the previous analysis Q to include perturbations both for the three- 
metric and for the scalar field associated with inflation. In particular, both the 
perturbations of the three metric and the scalar field (which represents matter 
and whose non-zero vacuum expectation value has been incorporated in the 
cosmological constant) are expanded in terms of scalar harmonics on the three 
sphere and just retained to second order []l(|. We then consider the resulting 
coupled equations in the presence of a positive cosmological constant (de Sitter) 
with the scope of obtaining information on the large-scale structures. 

Since we can not actually solve the nonlinear coupled matter-gravity system 
we do not know the effect of fluctuation corrections on the scalar field energy 
and the resulting distorted scalar field spectrum. However we can estimate 
the fluctuation corrections perturbatively using, as unperturbed solution, the 
matter wave function obtained on neglecting them (BO approximation). One 
then assumes that the distorted scalar field spectrum is reflected by the relative 
weight of the modes (unperturbed plus contributions due to fluctuation cor- 
rections) at a given instant, say when each mode exits the horizon during the 
inflationary era [JL2]. It is further clear that the relative strength of the diverse 
corrections besides depending on the moment at which they are frozen will also 
depend on the initial value of the cosmological constant, in any case the fluc- 
tuation corrections must always be less than the unperturbed result otherwise 
our perturbative approach is meaningless. 

In order to illustrate our approach we recall the perturbed minisuperspace 
results previously obtained |]Io| . The spatial metric h a p is given by 

h a (3 = a 2 (v) + Ea/s) j (1) 

where a is the Robertson-Walker scale factor, r\ the conformal time, VL a p the 
metric on a sphere of unit radius and e a p a perturbation to it which can be 
expanded in terms of the spherical harmonics of S 3 pi : 
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where the coefficients a n i m , b n i m are functions of 77 but not of the spatial coor- 
dinates. Similarly, one may also expand the lapse (N), shift (N a ) functions and 
the scalar (matter) field <f>: 
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and we have just considered terms in the expansions eqs. (||)-(||) associated 
with scalar harmonics Qf m on the three sphere (or their covariant derivatives 
[ p"l| ) since it is these contributions which are relevant for the non-homogeneous 
part of the scalar field in the semiclassical limit Jl(J . Henceforth, for the sake of 
brevity, we shall only exhibit the index n, other indices (I, m) being understood. 

One may now substitute the above into the total Hamiltonian density (grav- 
itation plus matter) and perform the spatial integration keeping the minisuper- 
space variable a to all orders and the perturbations </>, a n , b n , /„, k n only to 
second order. These now become the dynamical variables and one may now 
quantize canonically obtaining the following WD equation and momentum con- 
ditions: 



V n ) 



H u + H M ) * = 



"H^ty = 0. 



(6) 



(7) 



The indices 0, 1 and 2 indicate the order with respect to the perturbations, and 
S the scalar part of the total Hamiltonian (obtained by setting to the vector 
and tensor perturbations in E3| ) . H G and H M are respectively the gravitational 
and matter parts of the Hamiltonian and ^(a, {a ra }, {b n }, <f>, {f n }) is the total 
matter-gravity wave function. 

We now follow a procedure analogous to the one illustrated elsewhere M 
and factorize the wave function ^> as: 
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#(a , {a n } ,{&„}, 0, {/„}) = JJx„(a,a„ ,b n , 0, f n )i> n (a 

n 

= n^n^=^' ( § ) 

n I 

which on substituting into eq. ^ leads to the following coupled equations: 

(H G + (H M ))4> = -(Hg n )4>, (9) 

(H M - {H M ))x + ^(V log$ r GV* = -{Hg n - (HgM , (10) 

where HE is the gravitational kinetic energy term and the averaging procedure 
(} is over all matter field configurations. Further G is a matrix depending on 
the diverse gravitational degrees of freedom, V is a vector gradient with respect 
to the various gravitational degrees of freedom and V T its transpose. It is 
understood that eq. ( |Io|) is evaluated where \ip\ 2 has support 0. The terms on 
the r.h.s. of eqs. ([))) and ( [To|) are associated with fluctuations and disappear in 
the BO approximation. Analogous equations are obtained from the constraints 
eq- (0). 

We may now introduce a semiclassical approximation to the gravitational 
wave function through: 

xjj(a,{a n + b n }) = N G e* sG , (11) 

with 

iy a 2 -2Y S^bi , 1rA 

where S G is the classical gravitational action and is solution to the Hamilton- 
Jacobi equation obtained for h — > from eq. (^) in the absence of fluctuations, 
and the contribution to the lower order (in Ti) from Nq is negligible in our 
present approximation. Further, since the momentum constraints eq. ([?]) are 
associated with reparametrisation invariance and in general reduce the number 
of parameters leading to a mixing of matter and gravitation degrees of freedom 
| jl2"| in contrast with our factorization ansatz eq. (|b|), we require that they be 
satisfied just for the gravitational wave function ■ip in the absence of matter 
backreaction [Q. This then leads to the above parameter dependence in -0. 

Through the above, as previously explained @|-||, one may introduce a 
conformal time t] and eq. Jl(i| ) becomes: 
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H M ))x-in^x 



-(Hg n - {Hg n ))x (13) 



and in particular on neglecting the fluctuation in eq. (13) one obtains: 



which is the usual evolution equation for matter. Eq. (J14|) is then solved through 



an ansatz for Xs 13 



where: 



Xs = Ne* : 



(15) 



5" = S (a, 0)+ 



(16) 



which is substituted into eq. (|T^) . On equating coefficients of the same order in 
a>n ,b n , f n , a series of relations are then obtained and solved for the functions 

The solutions for ip and x are then used to obtain expressions of physical 
interest. Further, since a semiclassical limit for ip is considered, for the values 
of a, a n and b n one uses their average (classical) values. In particular, from the 
matter wave function Xs (solution to eq. ([l4])) one determines the expectation 
values of the coefficients f 2 at the exit from the horizon, obtaining with H the 
Hubble parameter: 



hH 2 

2n 3 ' 



(17) 



which is directly related to the spectrum of density fluctuations and corresponds 
to the Harrison-Zel'dovich (HZ) spectrum (since k = with k the physical 
wave number) py . Further, one may estimate, using the lowest order solution 
of eq. (|l4|), the magnitude of the fluctuations in eq. (HI]) through: 
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Figure 1: Unperturbed spectrum (HZ) and its distorsion (fluctuation corrections 
+/— ) power spectrum at ~ 10 17 GeV and ~ 10 14 GeV inflation energy. 



which may be interpreted as a correction to the scalar field energies and conse- 
quently as a distortion of the unperturbed (HZ) spectrum and we shall denote 
it by A(/^). One then obtains: 



= ±- 



(H kin^kin) (^kin) 2 



(19) 



{fl) (H M ) 

where the scale factor is evaluated at the exit from the horizon for each mode. 

In particular, if the spectrum is evaluated at a common time after having 
re-entered the horizon, one finds that the unperturbed part behaves as n, while 
the corrections behave as ~ n _1 and ~ n 2 for small and large n respectively and 
are therefore most effective both at extremely large and small scales. Further 
the interval in n for which our perturbative approach is valid increases as the 
cosmological constant driving inflation is decreased (see Fig. |l|). For example, 
for an inflation energy ~ 10 17 GeV and an e-folding of 100, the HZ spectrum is 
modified by our perturbation on physical scales in the interval (1-^10 4 ) h~ x Mpc 
(we adopt the value h = 0.5 for the present value of the Hubble constant Ho in 
units of 100 km s^Mpc^ 1 ). However, from the analysis of present inflationary 
models |H| one obtains that standard and chaotic inflation is consistent with the 
COBE results |§ g7| only for energies ~ 10 u GeV (which is the energy scale 
for which the observable fluctuations exit the horizon) and in order to obtain 
a higher energy one should construct models with suitable e-fol ding and/or a 
horizon growth during inflation steeper than that presently used |l5|. 

With the fluctuation spectra obtained (HZ ± fluctuation corrections) and 
for an inflation energy of ~ 10 14 GeV, we have performed numerical simulations 
assuming cold dark matter dominance (the so-called CDM model): the con- 
stituents of dark matter in this model arc massive particles, which decoupled 
from radiation when non-relativistic or never were in thermal equilibrium. 

In the last decade the standard CDM model has shown a high predictive 
power in explaining many observed properties of the large-scale galaxy dis- 
tribution. However, it is now known that this model has some serious prob- 
lems, mostly due to the high ratio of small to large-scale power. In particular, 
the COBE normalization jl6| |l7j implies excessive velocity dispersion on Mpc 
scales |l8| and is unable to reproduce the slope of the galaxy angular correlation 
function obtained from the APM survey fl9|| . The spectrum of the primordial 
fluctuations in the CDM model is the HZ spectrum and, as mentioned, our 
solutions modify this spectrum at small and very large scales. 

To follow the non-linear evolution after the matter-radiation decoupling 
we used a particle-mesh code ^0|. A preliminary analysis with N, p = 128 3 
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Table 1: Comparison of observed and computed mass variance. 



hr x Mpc 


APMgl) 


QF_ 


HZ 


10 


1.05-1.47 


1.10-2.08 


2.31 - 3.07 


15 


0.61 - 0.89 


0.58-0.85 


0.99 - 1.19 


20 


0.39-0.62 


0.35-0.43 


0.55-0.57 


25 


0.29-0.49 


0.21 - 0.24 


0.31 - 0.33 


30 


0.18-0.33 


0.13-0.16 


0.18-0.21 



particles and N g = 128 3 grid-points indicated that adding or subtracting our 
fluctuation corrections to the HZ spectrum in the large n region, increased or 
decreased respectively the mass excess or the bulk velocities with respect to 
the HZ results. Therefore, we proceeded with more detailed simulations with 
N p = 256 3 particles on N g = 256 3 grid-points, on a Cray T3D MCA 64-8. 

We ran three simulations, one with the HZ spectrum over a box of size 
128/i _1 Mpc and two with the spectrum obtained on subtracting the fluctuation 
corrections from the HZ spectrum (we shall denote it by QFS), one over a box 
of size \2i]h~ 1 Mpc and the other over a box of size hYlhT x Mpc. 

The amplitude of the primordial fluctuation spectrum is not determined 
by our free theoretical parameters (that is, the inflation energy and the e- 
folding number) but by the inverse of the rms mass fluctuation on a sharp- 
edged sphere of radius 8h~ 1 Mpc (as)- The COBE DMR detection of large 
angular scale anisotropies of the cosmic microwave background | fl6| then fixes 
the normalization and makes the model completely specified. 

A detailed study of the numerical results is under way, however a prelim- 
inary analysis of the variance of mass reported in Table [l] for the box of size 
128h~ 1 Mpc (the range of values for QF_ and HZ refers to the greatest and 
the least of three evaluations obtained following ^T|) gives an idea of the rele- 
vance of the fluctuation corrections to the development of large scale structure. 
Analysis of mock galaxy catalogues obtained by such simulations should then 
provide a quantitative answer concerning the possibility of solving the two major 
problems of the CDM model with quantum corrections. 

Further details of theoretical and numerical aspects will be presented else- 
where n ||. 
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